Міністерство освіти і науки України Національний університет „Львівська політехніка” ІКТА, кафедра «Захисту інформації» Звіт З КУРСУ “ КОМП’ЮТЕРНІ МЕТОДИ ДОСЛІДЖЕННЯ ІНФОРМАЦІЙНИХ ПРОЦЕСІВ ТА СИСТЕМ ” НА ТЕМУ: “ МЕТОДИ ЧИСЕЛЬНОГО РОЗВ`ЯЗУВАННЯ ДИФЕРЕНЦІАЛЬНИХ РІВНЯНЬ“ Варіант 2  Мета роботи – ознайомлення з методами чисельного інтегрування диференційних рівнянь. КОРОТКІ ТЕОРЕТИЧНІ ВІДОМОСТІ Диференціальним називається рівняння, в яке входять похідні невідомої функції. Приклад:

(1)
(2) Диференціальне рівняння (ДР), що містить лише одну незалежну змінну і похідні за нею, називають звичайними (ДР). Це, наприклад, рівняння (1). ДР, що містить декілька незалежних змінних і похідні за ними, називають рівняння в частинних похідних Порядком ДР називається найвищий порядок похідної (або диференціалу), який входить в рівняння. Звичайне ДР (ЗДР)
-го порядку в загальному випадку має незалежну змінну, невідому функцію та її похідні (або диференціал) до
-го порядку включно:
(3)
- незалежна змінна;
- невідома функція (залежна змінна);
- похідні цієї функції. Диференціальне рівняння
-го порядку, розв’язане відносно старшої похідної, може бути записано у вигляді:
(4) Щоб розв’язати ЗДР, необхідно мати значення залежної змінної та (або) її похідних при деяких значення незалежної змінної. Якщо ці значення задані при одному значенні незалежної змінної - така задача називається задачею з початковими умовами або задачею Коші. Якщо ці значення задаються при
або більше значеннях незалежної змінної - задача називається крайовою. Значення залежної змінної та її похідних називаються ще додатковими умовами, котрі в задачі Коші називаються початковими, а в крайовій задачі - граничними. Метод Рунге-Кутта-Мерсона з автоматичною зміною кроку Метод дозволяє оцінити похибку на кожному кроці інтегрування. Похибка інтегрування має порядок
. При цьому не потрібно зберігати в пам’яті обчислення значень функцій на кроці
і
для оцінки похибки - перевага порівняно з методом Рунге-Кутта з автоматичною зміною кроку. Алгоритм методу 1. Задаємо число рівнянь
, похибку
, початковий крок інтегрування
, початкові умови
. 2. За допомогою п’яти циклів з керуючою змінною
обчислюємо коефіцієнти

3. Знаходимо значення

та похибку
4. Перевіряємо виконання умов
Можливі випадки: а) Якщо перша умова не виконується, тобто
, то ділимо крок
на 2 та повторюємо обчислення з п.2, встановивши початкові значення
. б) Якщо виконується перша та друга умови, значення
та
виводяться на друк. Якщо друга умова не виконується, крок
збільшується вдвічі і тоді обчислення знову повторюється з п.2 (нема потреби обчислювати при малому кроці). Треба відмітити, що похибка
на кожному кроці методу Рунге-Кутта-Мерсона оцінюється приблизно. При розв’язуванні нелінійних ДР істинна похибка може відрізнятися в декілька разів від заданої
.
, де
.
- крок поділити на 2 і повернутися на початок.
для всіх рівнянь: виводимо на друк
, а крок збільшуємо удвічі. 2.ЗАВДАННЯ ДО ЛАБОРАТОРНОЇ РОБОТИ Розв’язати чисельним методом звичайне диференційне рівняння.

  СПИСОК ІДЕНТИФІКАТОРІВ КОНСТАНТ, ЗМІННИХ, ФУНКЦІЙ, ВИКОРИСТАНИХ У ПРОГРАМІ, ТА ЇХ ПОЯСНЕННЯ:                                                                   Текст програми на мові Паскаль: program runge_kut_merson; uses crt; const n=1;a=0;b=1;e=0.001; type vec=array[1..n] of real; var i,m,p :integer; x,h,k :real; y1,y,f,r,z :vec; Function dy(var k:integer;x,s:real):real; begin for i:=1 to n do z[i]:=y[i]+s; dy:=sin(sqrt(x...